The classical pion field in a nucleus. 
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Abstract 

A self-consistent symmetry arises when the nucleon angular momen- 
tum j and the isospin t are coupled to a grand spin G. Closed G shells 
become sources of a classical pion field with a hedgehog shape. Although 
the amplitude of the pion field, as measured by the chiral angle, is small, 
it is found to perturb significantly the energies of the nucleon orbits. 

1 Introduction. 

The purpose of this paper is to explore some nucleon configurations which give 
rise to a classical pion field. In most mean field calculations, the symmetries 
assumed for the intrinsic state do not sustain a classical pseudoscalar pion field. 
For example, a closed j shell can only give rise to a spherically symmetric mean 
field. The pseudoscalar classical pion field is expected to be composed of odd-/ 
spherical harmonics so that a closed j shell does not provide a source for the pion 
field. Nor do the usually assumed deformed intrinsic states of even-even nuclei. 
This is why only the exchange (Fock) term of the pion exchange interaction 
contributes to the mean field in these nuclei, while the direct (Hartree) term 
vanishes. 

This is a feature of mean field calculations in which the nucleus is described 
by a Slater determinant or a quasi-particle vacuum. However, it has been 
known for decades [1] that a considerable fraction of nuclear binding is due 
to 2-particle 2-hole excitations caused by the pion exchange interaction. Fur- 
thermore, Monte-Carlo calculations of light nuclei suggest that roughly 75% of 
the nuclear binding is due to the pion exchange interaction [2] . Mean field cal- 
culations, which use Skyrme or Gogny forces (see the review [3] and references 
therein) simulate the 2-particle 2-hole contribution of the pion exchange in- 
teraction by an effective density dependent interaction. Relativistic mean field 
calculations (see the review [4]) simulate the same effect by the presence of 
a sigma meson. However, the effective interactions (or mesons) used in such 
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mean field calculations are not able to explore the possibility that the pion field 
may favor somewhat different nucleon configurations which give rise to a finite 
classical pion field. 

Such a possibility has however been studied by Sugimoto, Ikeda, Toki et 
al.[5, 6, 7, 8]. They construct an independent particle state composed of either 
mixed parity nucleon orbits or mixed parity and isospin orbits. The charge and 
parity of this state are then projected and the energy of the projected state is 
minimized. They find that the tensor force, which is part of the pion exchange 
interaction, plays an important role, which is understandable insofar as the 
projection adds 2-particle 2-hole admixtures to the Hartree-Fock state. 

In this work, we investigate a self-consistent symmetry of the pion field which 
leads us to construct closed G-shells where G = j + t is the sum of the angular 
momentum and isospin of a nucleon orbit. 

2 The non relativist ic limit of Weinberg nucle- 
ons interacting with a chiral field. 

The lagrangian of Weinberg nucleons interacting with a chiral field U = exp {i^9 a T a ) 
is[9]: 

l = $ (id^ -m + i^m (fycr *)) i>+^tr [(d^u) (wrf)] +^y L tr (u + u*) 

(2.1) 

It can be written in terms of the chiral angle 9 a as follows: 

L = i> (id^ — M) ip 

+t (-e abc r a sin 2 B J± (l^) - 9a1s\ sin9 (7"^^) - 9aT*\^ (7^)) V> 

+\fl{[d,e){d^)+ S ^6[d^ (V^) -4m 2 sin 2 (2.2) 

where tp is the nucleon field, 9 2 — J2 a =i 3 an d where the last two terms in 
the second line have been multiplied by the axial coupling constant <ju[9]. The 
pion decay constant /„., the axial coupling constant gA, the pion and nucleon 
masses m w and M are respectively: 

U = 93.2 MeV g A = 1.257 = 139.6 MeV M = 938.9 MeV (2.3) 

Our results (see Fig. 3) yield chiral angles 9 a which do not exceed 0.05 so that we 
can safely work with the low amplitude limit, in which the chiral field reduces 
to a pion field. For small values of 9 a , the lagrangian (2.2) reduces to: 

L = i> [id^ -M- \ 9A i^T a {d^eS) i> + \ll (( W (&*o a ) - mlel) 

(2.4) 
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and the pion field 7r a is related to the chiral angle a by: 

7T a = UOa (2.5) 

The mean field approximation consists in assuming that the chiral (or pion) 
field is classical. We shall further consider time-independent stationary states, 
in which case we can work with the hamiltonian: 

H = Jd 3 x ^ (^-^- +pM- \g A ^r a (a ■ W0 a fj ^f 2 J d 3 x ((v0 a ) 2 + m 2 J 2 a 

(2-6) 

We further proceed to make a non relativistic reduction of the Dirac hamil- 
tonian: 



h 



a-V _ 1 



+ 0M- -9Al 5 r a (a ■ W a ) (2.7) 

The reduction is made in two steps. First we maintain the nucleon orbits or- 
thogonal to the vacuum Dirac sea orbits by expanding the eigenstates of h on 



the positive energy nucleon plane wave states k,a,r,+J, which are eigenstates 

of SL ¥- + pM belonging to the positive eigenvalues \Jk 2 + M 2 . This is achieved 
by working with a modified Dirac hamiltonian Iinr defined such that it has the 
same matrix elements as h between the positive energy plane wave states: 

k, a, r \h NR \ k', a', A = (k, a, r, + \h\ kf , a', r' , +) (2.8) 



The second step consists in expanding the resulting Iinr in powers of k/M. 
When the expansion is limited to order k 2 /M 2 the resulting non relativistic 
hamiltonian acquires the familiar form: 

hNR = ~-^<T i T a {V i B a ) (2.9) 

3 The model hamiltonian and G shells. 

In order to explore nucleon configurations which are likely to give rise to a clas- 
sical pion field, we further simplify the dynamics by assuming that the nucleons, 
which interact with the classical pion field, are in spherical oscillator orbits and 
that they also interact with a spin-orbit contact interaction of the form[3]: 

Vis=iWu X - J d 3 ri d 3 r 2 

V> f (ri) V f (f 2 ) \e ijk (a\ + cr l 2 ) (pj - p^ S (n - r 2 ) (p\ - p§)] V (r 2 ) V (n) 

(3.1) 

where p = V /i. The expectation value of the spin-orbit interaction Vi s in a 
Slater determinant |<I>) can be reduced to the form: 

($|V5a|$) - -4Wi s J d 3 r B (r) • (Vn (r)) (3.2) 
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where n (r) and B (r) are the following nucleon densities: 



n(f) = ($\^ (r)ip(r)\$) B (f) = ($ |V f (r) (a x p)^{r)\ $) (3.3) 
The model hamiltonian is: 

# = 1 <*W (h osc -^-a l r a (V i e a ))^ + V ls + ^ J d 3 x ((w a ) * + 

(3.4) 

where 



V 2 1 

2M 2 



(3.5) 



We shall diagonalize the single particle hamiltonian in a basis in which the 
angular momentum j and the isospin t are coupled to a grand spin G: 



\nljGM) = J2 (jm, \m t \GM) \nljm) 



1 

2 mt 



(3.6) 



where \nljm) is a usual spherical shell model orbit and (jm, ^m t \GM^ a 
Clebsch-Gordan coefficient. In this basis, the single particle spin-orbit inter- 
action remains diagonal. We will show that there is a self-consistent symmetry 
in which nucleons form closed G-shells and the pion field has a hedgehog shape. 
To any non-zero grand spin and projection (G, M) there correspond four states, 
which we denote as \GMi) with i = 1, 2, 3, 4, namely: 



\GM1) 


l=G-\=j-\ 


j = l + \ = G-\ 


P = {-l) u+l 


\GM2) 


I = G + 1 = j + \ 


j = l-\ = G+\ 


P = (-l) u+L 


\GM3) 


l = G = j-\ 


j = i + h = g + h 


P = (-lf 


\GM4) 


l = G = j + \ 


j = l-\ = G-± 


P = (-lf 



To grand spin G = states correspond two shell model states: 



|(G = 0),1) 


I = G = 


j — I 4" k — G ' + k — k 


P = + 


l(G = 0),2) 


Z=G+1=1 


j = l-\ = G+\ = \ 


P = - 



For example, the four \nGMi) states with n = 1 and G = 2 are: 



G 




Z 


j 




2 


1 


1 


3 
2 


1)D3 


2 


2 


3 


b 
2 




2 


3 


2 


b 
2 


Ids 


2 


4 


2 


3 

2 


Ids 



(3.7) 



(3.8) 



(3.9) 



and the two states with n = 1 and G = are the lsi and lpi states. 
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4 The self-consistent symmetry generated by closed 
G shells. 

A closed (nGi) shell is a state in which the 2G+ 1 states \nGMi) with different 
values of M are occupied. Applying Wick's theorem, we find that the expectation 
value of a tensor operator of rank K and component q in a closed (nGi) shell is 

(nGMi \T Kq \nGMi) = 5 Ka V2G+l (nGi \\T K \\nGi) (4.1) 

M 

so that it vanishes unless K — 0. When applying Wick's theorem, one should 
remember that the rank K of the tensor refers to rotations in the combined 
(r, cr, t) space which includes rotations in isospin space. The expectation value, 
in a closed G shell, of the classical pion field appearing in the hamiltonian (3.4) 
is: 

M 

To prevent £ from vanishing, r a a ■ V# a must be a tensor of rank zero in (far) 
space. Since r a is a tensor of rank 1 in this space, the condition for £ not to 
vanish is that a ■ V9 a should also be a tensor of rank 1. However, a ■ V is a 
scalar so that the condition requires 9 a to be a vector and this in turn means 
that the chiral angle must have the so-called hedgehog shape: 

a (r) = r a 6 (r) (4.3) 

where r — r/r. This also implies that closed G-shells have definite parity. 
If we calculate (Vi0 a ) with the hedgehog field (4.3), we find: 

^Ta {Vi? a 6) = a ■ t- + (f • 3) (f • r) ( ] (4.4) 

r \dr r J 

although straightforward calculation will also reveal that the two- 
orbit interaction V/ s gives rise to a single-particle potential of the 

h ls = m s -p-L + mJ- d ^- + ^) (4.5) 

r dr \r dr r A ) 

where L = fxp, and where n (r) and B (r) are related to the densities (3.3) as 
follows: 

n(r)=n (r) B (f) = ?B (r) (4.6) 
Each closed (nGi) shell makes the following contribution to these densities: 

n nGl (r) = ( nGMi 10 (r\nGMi) B nGl (r) = - (nGMi \f) -a-L (r\nGMi) 

M M r 

(4.7) 



nGMi J 



(4.2) 



A lengthy 
body spin 
form: 
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The total densities n (r) and B (r) are obtained by summing these contributions 
over the occupied (nGi) shells. Thus the single particle hamiltonian which gives 
rise to closed G shells is: 



hNR = + Im^v 1 - £± 

2M 2 2 



a • r — \- (r ■ a) (r ■ t) — 
r \dr 



+4W ls ~a ■ L + 4W ls 
r dr 



1 dB B 

r dr r 2 



(4.8) 



The expectation value of the third component t 3 = ^t 3 of isospin in a 
\nGMi) orbit is: 



M M 
(nGMi \t 3 \ nGMi) = — (i = 1, 4) (nGMi \t 3 \ nGMi) - ~ 2(G + 1) 



(i = 2,3) 
(4.9) 



It follows that a closed G shell has an equal average numbers of neutrons and 
protons. Furthermore, a closed G shell contains 2G + 1 orbits and it therefore 
contains an odd number of nucleons. 



5 Nucleon orbits in a hedgehog field of varying 
amplitude. 

A self-consistent calculation of the chiral angle requires an initial guess of either 
the nucleon orbits or of the chiral angle (r). In this respect, a useful guide con- 
sists in assuming a simple form for the chiral angle. Since it is a p wave, it likely 
to grow linearly with r at the origin. It is also likely to decrease exponentially 
at infinity as exp (—m^r). The simplest function which will do this is: 

9{x)=axe~ x x = m^r (5.1) 

where a measures the amplitude of the chiral angle. The self-consistent calcu- 
lation described in section 6 shows that (5.1) is quite a good initial guess for 
the shape of the chiral angle. An example is displayed in Fig. 3. We can then 
diagonalize the single particle hamiltonian (4.8) with and plot the spectrum of 
the nucleon orbits as a function of the amplitude a of the chiral angle (5.1). 

In this exploratory calculation, we do not mix radial quantum numbers so 
that each orbit has a given principal quantum number n. The nucleon orbits |A) 
are thus expanded on a \nGMi) basis: 

| A) = \n\G\M\k\) = \nxG x M x i) (5.2) 

i 

The sum over i extends from 1 to 4 when G > and from 1 to 2 when G = 0. 
However, since closed G shells give rise to orbits with definite parity only two 



G 



\nGMi) states contribute to any one orbit |A). The coefficients q 7 ^ x are 
determined by diagonalizing the hamiltonian (4.8): 

h N R \n x G x M x kx) = e nxGxkx \n\G\M\k\) (5.3) 

in this basis. The energy of filled G shells is then the following functional of the 
chiral angle 9: 

E(6)= ]T (2G A + 1) e nxGxkx + \fl J cPr (j^f + ^ + +($ \V ls \ $) 

(5.4) 

where the last term is the expectation value of the spin-orbit interaction (3.1). 
The sum is limited to occupied G shells. 

Figure 1 shows the spectrum of the low energy nucleon orbits plotted as a 
function of the amplitude a of the chiral angle. They are calculated using a har- 
monic oscillator constant fiui = 12.71 Mev which fits the observed charge radios 
2.71 fm of 16 O[10]. Figure 2 shows somewhat higher energy orbits calculated 
with huj = 10.27 MeV which fits the observed charge radius 3.48 fm of 4a Ca. 
In both figures the spin-orbit interaction strength was Wi s = 15 /m~ 5 adjusted 
to yield a 5 Mev splitting between the Ids and Idz states in 16 C When a = 
the chiral angle and the pion field vanish and the orbits reduce to spherical shell 
model orbits \nljm). In the presence of the chiral field (a ^ 0) each j shell of 
given parity splits into two G shells with G = j ± | . 



6 Self consistent calculation of the chiral angle. 



The chiral angle can be calculated by minimizing the energy (5.4). The equation 
SE (8) /59 (r) = can be solved for 9 (r) in the form: 

777 f°° 

6 ( X ) = -T^72 x ' 2dx ' ii(z<)M*>)p(*') (6-1) 

ZlT J-K JO 

where x < and x> are respectively the lower and upper values of x and x'and 
where the Bessel functions are: 

, , sinh (z) cosh (z) , . , tt „ / 1 \ . 
h (z) = 5^ + — fci {z) = —er z 1 + - (6.2) 



The source term in (6.1) is calculated in terms of the occupied orbits: 

x2p{x) = Wix) £ (2G x + l)e nxGxkx (6.3) 

The spectra displayed on Figs.l and 2 can be used to guess which filled G 
shells are likely to give rise to a favored configuration. We expect that favored 
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configurations will be those which give rise to an appreciable energy gap between 
the occupied and empty G-shells. Once we have decided which G shells are filled, 
we can make a self-consistent calculation of the chiral angle 9 (r) : 

1. We choose a set of closed G shells and, with this set, we calculate the 
source term p (x) using (6.3) with unperturbed orbits. 

2. With this p(x), we calculate 6 (x) using (6.1). 

3. With this 9 (x) we calculate the orbits (5.2) by diagonalizing the hamilto- 
nian (4.8). 

4. With this set of orbits we calculate p (x) using (6.3) and we return to step 
2. We continue iterating until convergence is reached. 

The results are summarized in Tables 1 and 2 for configurations containing 
12 < A < 37 nucleons and 2p — If shell configurations respectively. The tables 
specify which closed G shells are filled, the number of nucleons, the energy gap 
separating filled and empty orbits and the maximum value of the chiral angle 
9(x) for each configuration. An orbit labelled, for example, (ld 5 / 2 ) G = 2 
refers to an orbit which, on Fig. 1 , is marked G — 2 and which converges to the 
ld 5 / 2 state when a = 0. Although the self consistent orbits are not identical to 
those plotted on Fig.l, they have a large overlap with the latter. We have only 
listed cases where the filled G shells have lower energies than the empty shells. 

When both G shells emanating from a given j shell arc filled, the resulting 
Slater determinant is equivalent to one in which the j shell is closed and it does 
not contribute to the source term p (x) . The closed j shells are listed in the first 
column of tables 1 and 2. The second column specifies the filled and closed G 
shells which correspond to partially filled j shells. Since every G shell contains 
an odd number of nucleons, and since a closed G shell has the same average 
number of neutrons and protons, configurations which have an odd number of 
filled G shells refer to odd-even nuclei. Most configurations listed in table 1 are 
of this type. Thus, for example, the A = 13 configuration could refer to the 
mirror nuclei ? 3 N e and g 3 CV and the A = 31 configuration to fgSis and i5-Pi6- 

There is a difference between the energy gaps, displayed in Figs.l and 2, and 
the energy gaps which are given in Tables 1 and 2. The energy gaps appearing 
in the figures are plotted as a function of the strength of the pion hedgehog field 
(as measured by the amplitude of the chiral angle). However, in Tables 1 and 
2, the displayed energy gaps are obtained self-consistcntly by a minimization 
of the energy. We see that in most odd-even nuclei, the chiral angles which 
minimize the energy are considerably smaller than the ones which produce the 
largest energy gaps (displayed in Figs.l and 2). The reason that the energy 
required to produce a classical field increases with the amplitude of the chiral 
angle. This is the case of most odd-even nuclei, listed in Tables 1 and 2, for 
which the chiral field is not strong enough to produce a level crossing between 
filled and empty levels. 

One exception in Table 1 is the A = 22 configuration, which has filled 
(1^5/2) G — 2 and (2s 1 / 2 ) G = shells. Fig.l shows that, for a < 0, the 
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(2si/ 2 ) G = orbit crosses the (1^5/2) G = 3 orbit and it does so also with 
the calculated chiral angle. Although it has a low energy gap, this configuration 
gives rise to a relatively strong chiral field. It refers to an odd-odd nucleus such 
as llNan. 

Another example is the A = 62 configuration shown in Table 2. In this con- 
figuration, three filled G shells become almost degenerate, namely the (I/7/2) G = 
4, (2p 3 / 2 ) G — 2 and (2pi/ 2 ) G = orbits. This is clearly visible in Fig. 2. In 
the self-consistent calculation, they are also almost degenerate, their energies 
being respectively 53.33,54.68 and 54.19 MeV. This configuration also gives 
rise to a relatively strong chiral field. It refers to an odd-odd nucleus such as 
i?Go 3 i. 

The configurations containing A = 13, 15, 29, 31 and 37 nuclcons in Table 
1 and A = 62, 65 and 67 in Table 2 have appreciable energy gaps. For these 
configurations, the chiral angle reaches maximum values 0.03 < 6 max < 0.04. 
For the shape 6 (x) = axe~ x this corresponds to a = e9 max in the vicinity of 
±0.1. In fact, the shape axe~ x is quite similar to the calculated shape of 9 (x) 
as shown on Fig. 3 and the self-consistent single particle spectra are very similar 
to those which can be read off Figs.l and 2. 



7 Discussion and conclusion. 

In most nuclear structure calculations, the pi-nuclcon interaction is included 
in terms of pion exchange interactions between nucleons and the possibility of 
generating a classical pion field is usually neglected. In this calculation were have 
explored such a possibility by constructing independent particle states which 
generate a hedgehog shaped classical pion field. We have made a parameter free 
estimate of the generated pion field and we found it to be small: the amplitude 
of the pion field, as measured by the chiral angle, does not exceed 0.05. We have 
found however, that such a classical field can perturb significantly the spectrum 
of the nucleon orbits. Such effects cannot be detected when Skyrme and Gogny 
interactions are used or in most relativistic mean field calculations (see the 
reviews [3] and [4]). Indeed, in such calculations, the pi- nucleon interaction 
is absent and the effects of pion exchange between nucleons is only included 
in terms of effective interactions (often represented by scalar and vector fields) 
which simulate pion exchange in the Fock term and 2-particle 2-hole excitations. 
The present calculation shows that, in the mean field approximation, the pion- 
nucleon interaction does give rise to a classical pion field in even-odd and odd- 
odd nuclei and it points to possible deformations which could be taken into 
account in, for example, generator coordinate calculations. 

We have considered the simplest case of a hedgehog shaped pion field. We 
have shown that such a field gives rise to nucleon orbits which have grand spin 
G = J + T and that the hedgehog shape of the pion field is a self-consistent 
symmetry in nuclei composed of closed G shells (which have G = J + T = 0). 
Closed G shells are states with equal average numbers of neutrons and protons. 
A single closed G shell is therefore an even-odd system of nucleons and two 
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closed G shells form an odd-odd system. 

If such states were stable, they would generate J = T rotational bands, 
as they do for example in Skyrmions[ll] and chiral solitons (see the review 
[12]). However, because of the weak amplitude of the pion field, the energy gaps 
created by the pion field in open j shells are not very large and it is unlikely 
that such rotational bands exist. A rotating nucleus composed of closed G shells 
would be seriously distorted by Coriolis forces acting on the rotating intrinsic 
state [13]. 

The present calculation does not therefore allow a quantitative prediction of 
the effect of the classical pion field on the spectra of some odd-even and odd-odd 
nuclei. It is an exploratory calculation which indicates the kind of deformations 
that are likely to be important in, for example, generator coordinate calculations 
of the spectra of odd-even and odd-odd nuclei. We have limited the present 
calculation to closed G shells which generate a hedgehog shaped pion field. 
There is no reason why this self consistent symmetry should be favored in all 
nuclei. Other self-consistent symmetries may be energetically favored when, for 
example, the neighboring even-even nuclei have quadrupolc deformations. 



A Some properties of grand spin states. 

Simple expressions may be derived for matrix elements between grand spin 
states \nljGM) (3.6). If we adopt the conventions (3.7) and (3.8) to designate 
the four grand spin states \GMi) which correspond to a given (G, M), then the 
matrix elements of a ■ f are: 



/ 1 o 
1 



(GMi \a • f | G'M'j) = 5 gg ,8mm>M^ t) 
\ 











2G+3 
2G+1 
VG(G+1) 
2G+1 



VG(G+1) 
2G+1 

2G-1 
2G+1 



(G>0), M\ 



(<"-) 



The matrix elements of a ■ r are: 



{GMi \a • f| G'M'j) = 5 G g>6mm>M% 



(G>0), M^ a) 



(ra) 



( 











-1 












-1 










-1 










\ 


-1 











) 



-1 





The matrix elements of (<? ■ r) (r • r) are: 

(GMi | (a • f) (r • r)\ G'M'j) = &gg,6 mm ,M s 



(rarr) 



-3 
1 



(G = 0) 



(G = 0) 
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M (rcrrT) 



1 V(G(G+1)) 
2G+1 2G+1 
V(G(G+1)) 1 
Z 2G+1 2G+1 





2G+1 





M (r<TrT) 



1 o 

-1 



dWG(G + l) 
(G = 0) 



2G+1 






y/G(G + l) 
i 

2G+1 



The spin-orbit interaction L ■ s is diagonal: 



(GMi 



L-s 



- G'M'j ) = Sgg'Smm' Sij A 



i\{G-\). 



\{G + 2), 1 -G,-\{G + 1)) (G; 0) 



Af s) = (0,-1) (G = 0) 

It is sometimes useful to use another grand spin basis \nlDGM). It is 
defined by coupling the spin s and the isospin t to a total spin D = s + t and 
by coupling the orbital angular momenta Z and D so as to form the grand spin 
G = 1 + D. For a given (G, M) there are four such states when G > and two 
states when G = 0: 



|(D) GMI) 


I = G- 1 


D = 1 


p =( _ 1) G + l 


|(D) GM2) 


Z = G+1 


D = 1 


P= (-l) U + i 


|(D) GM3) 


Z = G 


D = 


P=(-lf 


|(D) GMA) 


Z = G 


D = 1 


P=(-lf 



(G>0) 



|(D) GMI) Z = D = P = +l 
|(D) GM2) 1 = 1 D = l P = -l 



(G = 0) 



(A.l) 



The grand spin states |(D) GMi) are related to the grand spin states \GMi) as 
follows: 

|GM1) = |(D) GMI) 
\GM2) = \{D)GM2) 



\GM3) 



G + l 
2G+1 



|(D) GM3) 



G 



2G+ 1 



|(D) GM4) 



|ZjGM4) = -yJ-JL-KD) GM3) + J |(D) GM4) 



(G>0) 



|G = 0,i) = |(D)G = 0,i) 



(A.2) 
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B The contribution of the spin-orbit interaction 
to the mean field. 

Using the v (1, 2) = v (2, 1) symmetry of the two body spin-orbit interaction, its 
second quantized form (3.1) becomes the sum of four terms: 

v ls = vw + y (2) + u (3) + u (4) (b.i) 

which are: 

U« = iW ls s ijk J d 3 raUl (a \a^\r) (r\p k \ S) ((3 \r) (r | 7 ) a 7 a s 

V ( V = iWi s s ljk J dld2al a\ (a \a l \r) (r\S) (/? (r\p k \ 7 ) a 7 a 5 

= -iW ls e ijk J dld2aial (a |aV| r) (r\S) (/? |r) (f \p k \ 7 ) a 7 a 5 

y( 4 ) =-^ s£ijfc y'dld2at e 4(a| ( T i |r7(r|/|5)(/3|^|r)(r|7}a 7 a 5 (B.2) 

In these expressions, |a) , |/3) , ... are nucleon single particle states expressed in 
an arbitrary basis and a sum is assumed over repeated indices. 

We evaluate the expectation value of VJ S in the Slater determinant: 

i*> = 4A---4j o > = IHio) (b.3) 

agf 

We use Wick's theorem to obtain: 

($ \V U \ $) = ($ |^i| $) + ($ \V 2 \ $) + ($ |V&| $) + ($ I V4 1 $) 



= 2iW, s 2 / d\ e ijk {(A |o-V| r) (r\p k \ A) (/x |r) (r + (A |a*| r) (r |A) (/x |^'| r) (f |/| /x) 

- (A |ctV| (^|A) (m |r> (r |p fc | it) - (A | cr* | f) (r\p k \ A) (/1 r) (r |/x}} 

(B.4) 

Using p k = iV fe , integrating by parts and using the antisymmetry of e^k we 
can reduce (<3? \ Vi s \ $) to the form: 

($ |V is | $) - 2W, a J] / d ' Ve ^ ( A 1*0 (r\^P k \ A) (V (a* |r) <r |^» 
+2W is J] / d 3 r £ ^( M |r)(r|p fc | A i)(V J (A|a 4 |r)(r|A)) 



-2IU is ^ y d 3 r £yfc (A |aV | r) {? |A) (/x |r> (r /i) 



A^G-F 
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-2W ls f &™Hk (A \r) (r |aV | A) (n \ip>\r) (f | M ) (B.5) 

This expression involves the densities: 

n{f) = Y, (A \r) (f | A) j? (r) = £ (A |f) (f |p] A) 

AGF AGF 

S(r) = ^(A|r)(r|(axp)|A) a(f) = £ (A |f) (f A) (B.6) 

AGF AGF 

in terms of which: 

($|V is |$) = 2VK Js y rf 3 r [-B(r) • (Vn (f)) -p(r) • (v x a(r)) - (r) • p (r) + iB ■ f (r) 

(B.7) 

Angular momentum coupling techniques can then readily be used to evaluate 
the contribution of a closed G shell to these densities, with the result that the 
densities have the form: 

1 - 1 ^.dn ^ ^1 

n(r) = n (r) p (f) = — Vn (f) = 7pr— 3 (f) = fir (r) B (r) = r-B (r) 

(B.8) 

which could have been be guessed from the fact that closed G-shells are invariant 
under rotations generated by spin + isospin. The second term of (B.7) therefore 
vanishes. The contribution of the spin orbit interaction to the energy of a closed 
G-shell nucleus reduces thus to: 



*|V5 a |$) = -4W ls J d 3 r B (r) • (Vn (r)) (B.9) 



( 

which is the expression (3.2). The contribution can be expressed in terms of the 
occupied orbits: 

($ \v u \ *) - -m s £ / d " r < A 10 x p)i\ A > ( v * & io W) (b.io) 

A^GF"' 

To obtain the contribution to the mean field, we calculate the (functional) 
derivative: 

" (^\V ls \^) = -4W l s[(rs\(axp) i \X)(V i n)-(fs\X)(V i B i )} (B.ll) 



5{\\rs) 

The contribution of the spin-orbit interaction to the mean field is thus the 
operator: 

his = -4W ls (V 4 n) (<? x p) i + W ls (V^) (B.12) 
Since n(f) — n(r), we have: 

1 dn _ . 1 dn _ ? 

Vin (cr x p) i = -— r ■ (a x p) = — a ■ L 

r dr r dr 
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V 4 B 4 = V^B =1^ + 4 (B.13) 
r ar r A 

so that: 

hu = m s -p • L + W ls (~ + ^) (B.14) 
r ar \r ar r l ) 

which is the expression (4.5). 



References 

[1] G.E.Brown. Unified theory of nuclear models and forces. North Holland, 
1967. 

[2] C. Pieper and R.B. Wiringa. Quantum monte carlo calculations of light 
nuclei. Annual Rev. Nucl. Part. Sci. 51, 2001. nucl-th/0103005. 

[3] Michael Bender, Paul-Henri Heenen and Paul-Gerhard Reinhard. Self- 
consistent mean-held models for nuclear structure. Rev.Mod.Phys. 75, page 
121, 2003. 

[4] D.Vretenar, A.V.Afanasjev, G.A.Lalazissis and P.Ring. Relativistic 
hartree-bogoliubov theory: Static and dynamic aspects of exotic nuclear 
structure. Phys. Rep. 409, page 101, 2005. 

[5] Hiroshi Toki, Satoru Sugimoto, and Kiyomi Ikeda. Relativistic mean-field 
theory with the pion in finite nuclei. Prog. Theor. Phys. 108, page 903, 2002. 

[6] Yoko Ogawa, Hiroshi Toki, Setsuo Tamenaga, Hong Shcn, Atsushi Hosaka, 
Satoru Sugimoto and Kiyomi Ikeda. Chiral sigma model with pion mean 
field in finite nuclei. Prog. Theor. Phys. Ill, page 75, 2004. nucl-th/0312042. 

[7] Yoko Ogawa, Hiroshi Toki, Setsuo Tamenaga, Satoru Sugimoto and Kiyomi 
Ikeda. Charge parity projected relativistic mean field model with pion for 
finite nuclei. Phys. Rev. C73, page 034301, 2006. 

[8] S. Sugimoto, K. Ikeda and H.Toki. Study of the effect of the tensor correla- 
tion in oxygen isotopes with the charge and parity projected hartree-fock 
method. Phys. Rev. C75, page 014317, 2007. nucl-th/0607045. 

[9] Steven Weinberg. The Quantum Theory of Fields. Cambridge University 
Press, 1986. 

[10] B.A.Brown, C.R.Bronk and P.E.Hodgson. Systematics of nuclear RMS 
charge radii. J.Phys.G:Nucl.Phys. 10, page 1683, 1984. 

[11] C.R.Nappi G.S.Adkins and E.Witten. Nucl.Phys. B228, page 552, 1983. 

[12] Georges Ripka. Quarks Bound by Chiral Fields. Oxford University Press, 
Oxford, 1997. 

[13] J.P.Blaizot and G.Ripka. Phys. Rev. D38, page 1556, 1988. 



14 




Figure 1: The spectrum of the \nGMk) orbits are plotted as a function of the 
amplitude a of the chiral angle (5.1). When a — the chiral angle vanishes and 
the orbits reduce to spherical shell model orbits \nljm). When a ^ each j 
shell splits into two G orbits with G = j ± | as indicated on the figure. The 
orbits were calculated with an oscillator frequency hui = 12.71 fm applicable 
to le O. 
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0.05 




-0.25 - 



Figure 3: The full line is the chiral angle 9 (x) and the dashed line is the source 
term p{x) defined in (6.3). They are plotted as a function of x = m^r for the 
A = 31 configuration. The pion field ir a (r) can be expressed in terms of 9 as 
n a (r) — f^r a 9 (r). In order to fit the figure, the source term p (x) is divided by 

47^/2. 
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huj = 12.71 MeV 










Closed j shells 


Filled G shells 




Gap (MeV) 


@max 


lsi/2 G = 0, 1 
1 P 3/ 2 G = 1,2 


(lp 1/2 ) G = 


13 


4.14 


0.038 


l*i /2 G = 0,1 
1 P 3/ 2 G = 1,2 


(lpi /2 ) G = 1 


15 


4.27 


-0.031 


lsi/2 G = 0, 1 
1 P 3/ 2 G = 1,2 
lp 1/2 G = 0, 1 


(l4/ 2 ) G = 2 


21 


1.17 


-0.017 


i*i/ 2 G = 0,l 

l?>3/2 G = 1,2 

lpi/2 G = 0, 1 


(l4 /2 ) G = 2 
(2*i/ 2 ) G = 


22 


1.4 


-0.04 


i*i/ 2 G = 0,l 
1 P 3/ 2 G = 1,2 
lpi/2 G = 0, 1 


(l4/ 2 ) G = 3 


23 


1.31 


0.02 


i*i/ 2 G = 0,l 

1P3/2 G = 1,2 

lpi/2 G = 0, 1 
ld 5 /2 <? = 2, 3 


(2*i /2 ) G = 


29 


2.51 


-0.035 


i*i/ 2 G = 0,l 
lps /2 G = l,2 
lpi/2 G = 0, 1 
14/2 G = 2,3 


(2*i /2 ) G = 1 


31 


2.92 


0.03 


hw = 10.27 MeV 










l*i /2 G = 0,1 
lp 3 / 2 G = l,2 
lpi/2 G = 0, 1 
14/2 G = 2, 3 
2si /2 G = 0,l 


(l4/ 2 ) G = 1 


35 


1.64 


0.26 


i*i/ 2 G = 0,l 
1 P 3/ 2 G = 1,2 
lpi/2 G = 0, 1 
14/2 G = 2, 3 
2*i /2 G = 0,1 


(l4/ 2 ) G = 2 


37 


2.55 


-0.023 



Table 1: Closed G shell configurations composed of 12 < A < 40 nucleons. The 
first column lists the filled and closed j shells which do not contribute to the 
source of the pion field. The second column lists the occupied filled G shells. 
The column marked A specifies the number of nucleons in the configuration. 
The last two columns list the energy gap separating filled and empty orbits and 
the maximum value of the chiral angle 8. 
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Huj = 10.27 MeV 










Closed j shells 


Closed G shells 


,4 


Gap (MeV) 


@max 


Is, lp, 2s, Id 


(I/7/2) G = 3 


47 


0.56 


-0.011 


Is, lp, 2s, Id 


(l/ 7/2 ) G = 4 


49 


0.53 


0.013 


Is, lp, 2s, Id 
l/ 7/2 G = 3,4 


(2ps /2 ) G = 1 


59 


1.58 


-0.029 


Is, lp, 2s, lc? 
l/ 7/2 G = 3,4 


(2 P 3/ 2 ) G = 2 


61 


1.15 


0.023 


Is, lp, 2s, Id 
l/ 7/2 G = 3,4 


(2^3/2) G = 2 
(2p 1/2 ) G = 


62 


1.71 


0.043 


Is, lp, 2s, Id 










l/ 7/2 G = 3,4 
2p 3/2 G = l,2 


(2p 1/2 ) G = 


65 


2.18 


0.025 


Is, lp, 2s, Id 
I/7/2 G = 3, 4 
2p 3/2 G = l,2 


(2p 1/2 ) G = 1 


67 


2.06 


-0.02 



Table 2: The same as table 1 for 2p— If shell nuclei. For all these configurations 
the Is, lp, 2s, Id shells are filled. 
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